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Abstract 

We test the plausibility that a Majorana fermion dark matter eandidate with a sealar mediator explains 
the gamma ray exeess from the Galaetie eenter. Assuming that the mediator eouples to all third generation 
fermions we ealeulate observables for dark matter abundanee and seattering on nuelei, gamma, positron, 
and anti-proton eosmie ray fluxes, radio emission from dark matter annihilation, and the effeet of dark 
matter annihilations on the CMB. After disearding the eontroversial radio observation the rest of the data 
prefers a dark matter (mediator) mass in the 10-100 (3-1000) GeV region and weakly eorrelated eouplings 
to bottom quarks and tau leptons with values of 10“^-1 at the 68% eredibility level. 
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I. INTRODUCTION 


Since 2009 an increasingly significant deviation from background expectations has been iden¬ 
tified in the data of the Large Area Telescope (LAT) on board the Fermi Gamma Ray Space Tele¬ 
scope satellite [1-9]. The deviation appears around 2 GeV in the energy spectrum of gamma ray 
flux originating from an extended region centered in the Galactic Center. The source of the excess 
photons is unknown. Their origin can be dark matter (DM) annihilation, a population of millisec¬ 
ond pulsars or supernova remnants [10-13], or cosmic rays injected in a burst-like or continuous 
event at the galactic center [14]. It is, however, challenging to explain the excess with millisecond 
pulsars [15, 16] based on their luminosity function. 

Recently, several groups including Daylan et al. [7], Galore et al. [8], and the Fermi Collab¬ 
oration [9] re-analyzed data from the Fermi-LAT [17] and concluded that the 1-3 GeV gamma 
ray signal is statistically significant and appears to originate from dark matter particles annihi¬ 
lating rather than standard astrophysical sources. The peak in the energy distribution is broadly 
consistent with gamma rays originating from self-annihilation of dark matter particles [7, 18-23]. 
The intensity of the signal suggests a dark matter annihilation cross section at thermal freeze out 
[24-29]. The diffuse nature and morphology of the gamma ray excess is consistent with a Navarro- 
Frenk-White-like Galactic distribution of dark matter [8]. This gamma ray excess thus drew the 
attention of a number of particle model builders and phenomenologists [10, 14, 24, 30-32]. 

The conclusion that we have discovered dark matter particles, however, cannot be drawn yet. 
First, we have to be able to exclude the possibility of a standard astrophysical explanation. Second, 
we need to demonstrate that a dark matter particle that explains the gamma ray excess (with a given 
mass, spin, and interaction strength to the standard sector) is consistent with a large number of 
other observations. The latter concerns our paper. We aim to determine the microscopic properties 
of the dark matter particle from the gamma ray excess and check that these properties comply 
with limits from other experiments. We use dark matter abundance and direct detection data, 
measurements of the gamma ray flux from the Galactic Center, near Earth positron and anti-proton 
flux data. Cosmic Microwave Background (CMB) observations, and measurements of galactic 
radio emission as experimental constraints. 

Amongst the above listed experimental bounds the constraining role of radio emission has been 
debated in the literature. Bringmann et al. have shown that radio emission imposes severe con¬ 
straints on dark matter annihilation in the Galactic Center [33]. Radio emission, however, could 
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be induced by various processes including synchrotron radiation, inverse Compton scattering, ion¬ 
ization, and bremsstrahlung. Most studies of the radio constraint on dark matter, including that 
of Bringmann et ah, ignore energy loss processes other than synchrotron radiation. However, as 
pointed out by Cholis et al. in Ref. [34], there are several reasons why the other processes could be 
important. Cholis et al. have shown that after considering inverse Compton scattering induced by 
high densities of radiation in the inner Milky Way the radio constraint on dark matter is weakened 
by about three orders of magnitude [34] . As a result dark matter annihilating at the thermal rate 
remains compatible with the radio data. After considering the effect of diffusion the constraint 
will be further weakened. Due to this, we will exclude the radio data point from our combined fit. 

As theoretical description of dark matter we use the simplified model framework. Within this 
ansatz we make minimal and general theoretical assumptions. We consider a single dark matter 
particle that couples to various standard fermions via a mediator. Our dark matter particle thus 
annihilates to several final states which all contribute to the observables mentioned above. 

This paper is organized as follows. In Sec. II we introduce the simplified dark matter model we 
use. In Sec. Ill, we describe the observables of dark matter abundance and scattering on nuclei, 
gamma, e+, and p cosmic ray fluxes, and the effect of dark matter annihilations on the CMB. 
Our numerical results are given in Sec. IV. Finally in Sec. V we summarize our main results. We 
collect the formulae of Bayesian inference and likelihood functions in the Appendix. 


II. THEORETICAL HYPOTHESIS 

In this section we motivate and describe the theoretical hypothesis we test. In Ref. [35] we 
compared Bayesian evidences for three leading simplified models to explain the gamma ray ex¬ 
cess from the Galactic center. We found that the experimental data, especially dark matter direct 
detection, clearly preferred a Majorana fermion dark matter particle coupled to Standard Model 
(SM) fermions via a real scalar. Motivated by this we assume that the dark matter particle is a 
Majorana fermion, which we denote with x- Inspired by the Higgs portal mechanism [36], we use 
a simplified model to describe interactions between x and SM matter. We assume that the dark- 
standard mediator is a real scalar field, S, and the form of the dark matter to mediator coupling 
is 

D ( 1 ) 


3 


The presence of 75 is essential since it is lifting the velocity suppression that one otherwise en¬ 
counters in the indirect detection cross section, thus making this operator capable of explaining 
the gamma ray excess. The interaction between the mediator and SM fermions / is assumed to be 

Cs^XflfS. (2) 

In line with minimal flavor violation [37], we only consider the third generation fermions, i.e. 

/ = 

For simplicity we assume that mediator pair final states are not present in the dark matter 
annihilation and only consider s-channel annihilation diagrams. According to power counting of 
the dark matter transfer momentum or velocity [38], with the bi-linears in Eqs. (1) and (2) the 
annihilation cross section of the fermionic dark matter candidate is not velocity suppressed, that 
is cm ~ 1. The dark matter-nucleon elastic scattering cross section is spin-independent (SI) and 
momentum suppressed. 


III. OBSERVABLES 

In this section we describe the calculation of the observables that we use to constrain the pa¬ 
rameter space of our hypothesis. TABLE I summarizes these observables. 

A. Dark matter abundance 

We assume that dark matter particles, as standard thermal relics, have frozen out in the early 
universe acquiring their present abundance. We calculate this abundance using micrOmegas ver¬ 
sion 3.6.9 [45]. We imagine that x is the only dark matter candidate, that is we use a Gaussian 
likelihood function with a mean and width determined by PLANCK [39] 

= 0.1199 ± 0.0027. (3) 

It is challenging to estimate the theoretical uncertainty of the abundance calculation in a simplified 
model and the task is the subject of a separate paper. In supersymmetric models, for example, the 
theoretical uncertainty is comparable to the experimental one over the bulk of the parameter space. 
Based on this, we assume an extra theoretical uncertainty of the same size as the experimental 
error. 
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observable name 

expression 

experiment 

data points 

data source 

dark matter abundance 


PLANCK 

1 

Ref. [39] 

7 -ray flux 

dEdn 

Fermi-LAT 

24 

Ref. [8] 

cosmic e'''-ray flux 

d^e+ 

dE 

AMS-02 

72 

Ref. [40] 

cosmic p-ray flux 

d^p 

dE 

PAMELA 

23 

Ref. [41] 

Cosmic Microwave Background 

feS 

PLANCK 

3 

Ref. [42] 

dark matter direct detection 

S 

LUX 

1 

Ref. [43] 

radio emission 

Su 

Jodrell Bank 

1 

Ref. [44] 


TABLE I: Summary of observables we use to constrain our dark matter scenario. The expressions in the 
second column are defined in the text of this section. 

B. Gamma ray flux from the Galactic center 


In the theoretical scenario under scrutiny the excess gamma ray flux observed by Fermi-LAT 
is generated by the self-annihilation of x particles. The differential flux of photons as the function 
of energy E and observation region is given by 


dEdVt 




■m E Bs 

f 


dNf 
dE ■ 


( 4 ) 


Here {av) is the velocity averaged dark matter annihilation cross section at the Galactic center, 
Bf = {av)f/{av) is the annihilation fraction into the // final state, and dNj^/dE is the energy 
distribution of photons produced in the annihilation channel with final state //. The J factor in 
Eq. (4) is a function of the direction of observation -0 


J{^)= [ Pl{r)dl, 

J los 

with 


( 5 ) 


r = 
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The dark matter distribution in the Galaxy is deseribed by a generalized Navarro-Frenk-White 
(NFW) dark matter profile [46] 


Px(^) 


Po 


(1 + r/rsY~^ 


(7) 


Here = 20 kpe is the radius of the galaetie diffusion disk, Tq = 8.5 kpc is the so¬ 
lar distance from the Galactic center, and po is set to reproduce the local dark matter density 
Pxi'f'o) = 0.3 GeV/cm^. Following Refs. [7, 8] we fix the inner slope of the NFW halo profile to 
7 = 1.26 and set -0 = 5° in order to avoid bremsstrahlung and other secondary processes [24]. 

The differential yield dN^ /dE is different for the three final states we consider. As seen from 
Eq. (4), the total differential yield determining the gamma ray flux is the annihilation-fraction- 
weighted sum of the differential yields into specific final states. We sum over the contributions of 
the three individual SM fermions (6, t, r). As Bj depends on the model parameters, the gamma 
ray data plays an important role in constraining the coupling of the mediator to SM fermions. 

We use micrOmegas version 3.6.9 to evaluate the theoretical prediction for the differential 
gamma ray flux [45]. The gamma ray spectral data points that we input into our Gaussian likeli¬ 
hood function are taken from Ref. [8], including both statistical errors and empirical model sys- 
tematics. 


C. Cosmic positron flux near Earth 

The third generation fermion states produced by dark matter annihilation in our model can 
produce stable leptons in a variety of ways, including production via the decay of top quarks or 
tau leptons, or secondary production from hadron decays. These charged particles provide extra 
sources of cosmic flux in addition to the expected astrophysical backgrounds. Consequently the 
measurement of the electron and positron flux allows us to set constraints on the dark matter prop¬ 
erties. The predicted extra electron and positron flux from dark matter annihilation is comparable 
to that of the secondary production of electrons and positrons, which is one order of magnitude 
smaller than the measured electron flux itself. Thus dark matter annihilation affects more the 
positron flux than the electron flux. Since the prediction of the electron flux poses an additional 
challenge and it is the source of considerable uncertainties, we only focus on the positron flux and 
do not consider the electron flux or the positron to electron fraction in this paper. 

The propagation of positrons within the Galaxy is well-described by the following simplified 
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transport equation 


( 8 ) 


dfe+ 

dt 


V{K{E,r)VU) 


A 


(6(E,r)/e+) = Qe+{E,r), 


in the diffusion zone approximated by a cylinder with thickness 2L. In the above equation 
/e+(r, t, E) is the number density of positrons, K{E^ r) is the diffusion coefficient which is pa¬ 
rameterized as K{E, r) = Kq{E/GEVY, and b{E, r) is the rate of energy loss. The source term 
reads as 


Qe+(^,r) 


2m^ ^ dE j ' 


(9) 


with (iiVg+ /dE being the energy spectrum of positrons produced in the annihilation channel into 
//. The differential positron flux is given by 


d^e+ 

dE 



( 10 ) 


with Ve+ being the positron velocity. 

For the dark matter induced positron flux calculation in micrOmegas, we take the MED model 
for the above diffusion parameters: the index of the diffusion coefficient 5 = 0.7, the normalization 
factor Kq = 0.0112 kpc^/Myr, and the thickness of the diffusive cylinder L = 4 kpc [47]. For 
the astrophysical backgrounds, we adopt the following parametrization for the interstellar positron 
flux and the flux at the top of the atmosphere (TOA) [48] 


= C,+ E-^^+ + CgE-^^exp{-E/Eg), 

( 11 ) 


( 12 ) 


with best-fit parameters Ce+ = 72 s~^ sr“^ m“^ GeV“^, 7 e+ = 3.7, Cg = 
1.6 sr“^ m“^ GeV“^, 7 ^ = 2.51, Eg = 1 TeV, and solar modulation parameter 0e+ = 

0.93 GV obtained in Ref. [48]. In the above parametrization the potential obtains two contri¬ 
butions, one from the collisions of cosmic rays in the interstellar medium and another from the 
interactions of high-energy photons. 

As experimental input for the positron flux we use the new release of AMS-02 data [40]. We 
assume the theoretical uncertainty is the same as the AMS-02 experimental error and the form of 
the likelihood is a composite Gaussian [49, 50]. 
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D. Cosmic anti-proton flux near Earth 


The propagation of anti-protons through the Galaetie eylinder follows a similar diffusion equa¬ 
tion as Eq. (8) but there is an additional effeet from the galaetie wind and the souree term ineludes 
the annihilation of anti-protons in the interstellar medium as well as the annihilation of dark matter. 
The energy loss of anti-protons, however, is negligible compared with that of the positrons. 

The astrophysical background is calculated by adopting the set of propagation parameters 
called the KRA model in Ref. [51]. These parameters were extracted from a fit to standard cosmic 
ray data. To calculate the anti-proton flux from dark matter annihilation, i.e. d^p/dE, we use the 
same MED model described in the last subsection and assume the velocity of the convective wind 
to be Econv = 12 km/s. The anti-proton flux at low energies is also altered by solar modulation 
effects. We use the Eisk potential 0^?, which relates the local interstellar anti-proton flux to the one 
measured at the top of the atmosphere, as described in the KRA model. 

We use the latest release of PAMEEA data as experimental input for the anti-proton flux [41]. 
Note that the error bars in this data release are only statistical. Systematic error bars are expected 
to be of the same order of magnitude as in the first release of PAMEEA data [52]. We com¬ 
bine the uncertainties in quadrature and assume that the theoretical uncertainty is the same as the 
experimental error in the composite Gaussian likelihood function. 


E. Cosmic Microwave Background 


Dark matter annihilation in the early universe affects the CMB temperature and polarization 
fluctuations. Thus the CMB power spectrum measurement from PEANCK provides constraints 
on dark matter properties. A key quantity for determining the constraint on a given dark matter 
model is the efficiency for producing ionizing radiation. The authors of Ref. [42] provide values 
of the effective efficiency fes for different annihilation channels and dark matter masses that can 
be easily interpolated. We quantify the CMB constraints using the following likelihood function 


= exp 


-ip A 

2^"® ^ ^2 X 10-27cm3s-i 


GeV 




feS ^ ^ feS,i^ij 

i=T,b,t 


(13) 

(14) 


with Ai = 3.16 and ci = 4.64 for the PEANCK data. Here Bi is the annihilation branching 
fraction defined earlier. 
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F. Dark matter direct detection 


Direct detection of dark matter is facilitated by dark matter particles scattering on nuclei of a 
target material in a well shielded detector. The differential recoil rate of dark matter on nuclei, as 
a function of the recoil energy Er, is 

dR 


dE 


R m^rriA 


dvvf(v)^. 


R 


where is the nucleus mass, f{v) is the dark matter velocity distribution function and 

^ + (/I - 

Cl-^R -^n 


^max 


TT 


(15) 


(16) 


with Emax = 2p\v^/mA, G^{q^) = ^ [30], and = ^gsNN {N = p,n). EA{q) is the 

nucleus form factor and pA = (^x + is the reduced dark matter-nucleon mass. We 

assume that the local disk rotation speed is 220 km/s with the same value for the most probable 
speed of the dark matter’s Maxwell-Boltzmann velocity distribution. The Galactic escape speed is 
544 km/s [43]. 

As we only consider the interaction mediated by the scalar between the dark matter particles 
and the third generation quarks, the strength of the mediator-nucleon (N) interaction reads 


9snn — ^ruNfrc —• 

27 mf 

f=b,t ^ 


(17) 


Above fro = I - /r„ - /£ - fn and we adopt /£ = = 0.02, /£ = = 0.026, 

/t, = 0.043 [30, 53, 54]. 


For the LUX likelihood function, we use a Poisson distribution in the observed number of 
events N, 

(1) + s)N f,-ib+s) 

C{s\N) = P{N\s) = ^ -, (18) 

where b is the expected number of background events, 

roo J TD 

s = MT dE4>{E)—{E) (19) 

Jo “A/j 

is the expected signal, MT is the detector mass x time exposure, and (j){E) is a global efficiency 
factor that takes into account trigger efficiencies, energy resolution, and analysis cuts. Likelihood 
calculations are performed using a version of LUXCalc [55] modified to include the additional 
momentum dependence in Eqn. (16). For the LUX analysis region used by LUXCalc, N = 1 and 
b = 0.64; the efficiency curve (j){E) was generated by TPCMC [56] using the NEST model [57, 58]. 
See Ref. [55] for further details. 
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G. Radio Emission 


Electrons and positrons from dark matter annihilation are expected to lose energy through syn¬ 
chrotron radiation in the presence of large scale magnetic fields. Thus the radio emission in galax¬ 
ies and galaxy clusters can also be used to place constraints on the dark matter properties. The 
synchrotron flux density is given by 




4:71 pI 


dE du 


( 20 ) 


where d^f.+ /dE is the positron flux in units of (GeV cm^ s sr) The synchrotron power per 
frequency reads 


dWsyn \/3 e^B 


du bvr rrip 




■'syn 


( 21 ) 


with 


Stt 


The (5-function implies 


F(x) = xj K,„{0d( ~ ^S(x - 1 / 3 ), 


3eBp‘^ ^ 3eBE^ 
dvrm^ 47rm^ 


( 22 ) 


(23) 


For simplicity, we fix the magnetic field strength at a conservative lower limit [33] 


B = 50 pG. 


(24) 


The integration cone in the J factor corresponds to a 4" region around the Galactic center. 


IV. RESULTS 

We coded the Lagrangian of the relevant simplified dark matter model in FeynRules [59]. Cal¬ 
culation of observables, including the dark matter relic density and nucleon scattering interactions, 
differential gamma ray, and p fluxes, and radio signal were performed using a modified version 
of micrOmegas 3.6.9 [45]. Nested sampling and posterior distribution calculations were performed 
by MultiNest [60]. The nested sampling algorithm was developed to calculate marginalized pos¬ 
terior probability distributions and it is a Bayesian’s way to numerically implement Febesgue 
integration [61]. Since the relevant part of the likelihood distribution spans multiple orders of 
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parameter 


ms 

Aft 

Xt 

A. 

(unit) 

(GeV) 

(GeV) 




scan range 

1 - 10^ 

1 - 10^ 

10"^ - 10 

10"® - 10 

10"® - 10 

prior type 

log 

log 

log 

log 

log 


TABLE II: Scan ranges and prior types used for the scanned parameters. 



loglo(/n;^./GeV) 



FIG. I: Posterior probability distributions marginalized to the scanned model parameters. The likelihood 
function for these plots only contains the dark matter abundance and the anomalous Fermi-LAT gamma ray 
data. The dark and light regions hereinafter correspond to 68% and 95% credible regions, respectively. 

magnitude, we use log priors for all parameters. We present further details of our statistical anal¬ 
ysis in the Appendix. 

In our numerical calculation we fix the dark matter to mediator coupling as = 1, and we 
scan the following free parameters: 

P = (25) 

Here is the mass of the dark matter particle, ms is the mass of the scalar mediator, and A/ (/ = 
r, b, t) is the coupling of the mediator to the SM fermion pair ff as defined in Eqs. (1) and (2). 
The range of our scan over the above parameters and the type of prior we use is given in TABLE II. 
To build some intuition, first we examine the constraining effect of each observable one by 
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one. To this end we plot the posterior probability distributions marginalized to the seanned model 
parameters sueh that the likelihood funetion only eontains the dark matter abundanee and one of 
the other observables. In Fig. 1 we show marginalized posterior probability distributions taking 
into aeeount the dark matter abundance and the gamma ray data. The first frame of Fig. 1 confirms 
that the gamma ray data restrict the range of the dark matter mass close to 35-60 GeV [7, 8]. It is 
less appreciated, however, that uncertainties still allow a 40-100 (25-160) GeV dark matter mass 
range at the 68% (95%) credibility level. The gamma ray data, coupled with the dark matter relic 
density, allows the whole mass range of the scalar mediator. It is also interesting to note that the 
preferred dark matter mass region is dissected by a diagonal band with a lower posterior around 
the on-shell resonance region ms = 2m^. In this valley, dark matter resonantly annihilates via 
the s-channel mediator, depleting its abundance. Thus, it is harder for the model to match the 
PLANCK constraint. 

The relevant interaction strengths also remain virtually unconstrained as shown by the right 
frame of Fig. 1. PLANCK and the anomalous Fermi-LAT gamma ray data only restrict these 
coupling in the A., = 1 x lO-*^ - 2.5 x (6 x - 1) and Xt, = 2.5 x - 2.5 x 
(1.6 X 10“^ - 1) ranges at the 68% (95%) credibility level. Simultaneous order 1 couplings are 
marginally allowed and appear in the part of the parameter space where the annihilation cross 
section is suppressed by sizable ms- 

In Fig. 2 we show marginalized posterior probability distributions with the likelihood function 
containing only the dark matter abundance and the AMS-02 positron flux data. The AMS-02 
measurement of the positron flux features a small upward kink, a sudden change of slope, around 
35 GeV. A smooth background prediction has a hard time to reproduce this kink and systematically 
falls below the experimental points in the 35-50 GeV region. Positrons originating from the 
annihilation of a 35-50 GeV dark matter particle can fill the gap between the background and the 
data. Hence the AMS-02 data show a mild preference toward a dark matter candidate with 16-65 
(10-160) GeV mass at 68% (95%) credibility level. AMS-02 also restricts the dominant decay to 
r leptons with a A,- coupling around 2.5 x 10“^ — 0.1 and A;, coupling below about 10“^ at the 
68 % credibility level. 

Fig. 3 shows marginalized posterior probability distributions with the likelihood function in¬ 
cluding the PLANCK and PAMELA anti-proton cosmic ray flux data only. The PAMELA data in 
itself does not prefer any particular parameter region. Dark matter and mediator masses are both 
allowed in the full scanned range at the 95% credibility level. This happens with the exception 
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logio(;M;^,/GeV) 


log,o(Afc) 


FIG. 2: Posterior probability distributions marginalized to the scanned model parameters. The likelihood 
function for these plots only contains the dark matter abundance and the AMS-02 positron flux data. 

of a small island around ~ ms ~ 10 GeV where the combined PLANCK and PAMELA 
constraints are harder to satisfy. The reasons for this are that this island falls on the ms = 2m^ 
resonant annihilation corridor and the PAMELA data around 10 GeV leave very little room for 
dark matter. This situation improves for lower dark matter masses. 

Since dark matter masses above the top quark mass are allowed by the combination of 
PLANCK and PAMELA, the Xt coupling comes into play. These data, however, are not suffi¬ 
cient to constrain A^. It is interesting to note that PLANCK and PAMELA allow fairly large values 
of At, Aft, and A^ for heavier m^ and ms in order to accommodate the correct relic abundance. 

In Eig. 4 we show marginalized posterior distributions for the dark matter abundance and CMB 
likelihood function, defined in Eq. (14). The diagonal depletion of the likelihood function due 
to dark matter resonant annihilation is apparent in the m^ vs ms frame. The CMB l ik elihood 
function suppresses the posterior around m^ = 10 GeV providing more constraint on low mass 
dark matter. 

The posterior probability distribution projected to the Xb vs. A,- couplings shows a peculiar 
pattern. This pattern is the combined result of two relatively simple sets of constraints. Dark 
matter abundance is responsible for the low likelihood values at low Aft and A^. It turns out that 
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log,o(;«;^,/GeV) 
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FIG. 3: Posterior probability distributions marginalized to the scanned model parameters. The likelihood 
function for these plots only contains the dark matter abundance and the PAMELA anti-proton cosmic ray 
flux data. 

the PLANCK constraint on the amount of relic dark matter is hard to respect unless one of these 
couplings is sizable, that is Xb or A,- > 10“^ at the 68% credibility level. If both of these couplings 
are small then annihilation is slow and dark matter is overproduced in the early universe. In the 
large coupling region, on the other hand, the CMB constrains Xb and A,- from above. If any of 
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logioCAfo) logio(Afo) 


FIG. 4: Posterior probability distributions marginalized to the scanned model parameters. The likelihood 
function for these plots only contains the dark matter abundance and CMB terms, as defined in Eq. (14). 

these eouplings are larger than about 0.1 then dark matter tends to beeome under-produeed and 
the CMB reeeives too mueh modifieation from dark matter. The Xt eoupling is hardly eonstrained 
by the CMB at the 95% credibility level. 

In Fig. 5 we show the marginalized posterior distribution with the likelihood function contain¬ 
ing only the PLANCK and LUX data. Due to momentum suppression of the nucleon-x elastic 
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FIG. 5: Posterior probability distributions marginalized to the scanned model parameters. The likelihood 
function for these plots only contains the dark matter abundance and LUX data. 

scattering, the LUX data very weakly eonstrain the dark matter or mediator mass. In the vs. 
mg plane the diagonal resonant annihilation valley is visible, but no other structure is present. The 
posterior probability distribution for the couplings is very similar to that in Fig. 4. Similarly to the 
ease of the CMB, PLANCK and LUX only impose a eonstraint on the order 1 eouplings. 

As discussed in the Introduetion, the radio signal potentially very strictly constrains dark mat- 
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FIG. 6: Posterior probability distributions marginalized to the scanned model parameters. The likeli¬ 
hood function for these plots contains all of dark matter abundance, Fermi-LAT gamma ray data, AMS-02 
positron flux data, PAMELA anti-proton ray flux data, CMB and LUX data. 

ter [33] . Assuming that dark matter eontributes to the radio signal only by synehrotron radiation we 
find the radio flux upper limit of Jodrell Bank at 408 MHz [44] exeludes the dark matter hypothesis 
we eonsider by two orders of magnitude. Our finding fully eonfirms that of Ref. [33]. This exelu¬ 
sion, on the other hand, is lifted if inverse Compton seattering, ionization, and bremsstrahlung are 
also eonsidered as potential dark matter energy loss mechanisms leading to radio emission [34] . As 
shown by Ref. [34] the bound from the radio data is weakened by about three orders of magnitude 
if inverse Compton scattering is considered and is expected to pose no constraint after including 
Galactic diffusion effects. Due to this, we do not include the radio emission data point in our 
combined likelihood function. 

The summary of all constraints is presented in Fig. 6. The constraint on the dark matter mass 
is dominated by the gamma ray data and the final combination restricts to the 10-100 (7-125) 
GeV region with 68% (95%) credibility. Less of the low mediator mass region survives the scrutiny 
of the combined constraints, leaving the 3-1000 GeV ms region preferred at the 68% credibility 
level. The combined constraints prefer a somewhat correlated pair of A,- and A;, couplings in the 
intermediate 10“^-1 region at the 68% credibility level. Small (A < 10“^) and large (A > 1) 
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values of couplings are disfavored mostly by PLANCK at the 95% credibility. Comparing the 
second frame of Fig. 6 to those showing the individual constraints it is clear that simultaneous 
order 1 couplings are mildly under stress from almost all the data. 


V. CONCLUSIONS 

In this work we perform a comprehensive statistical analysis of the gamma ray excess from the 
Galactic Center in a simplified dark matter model framework. According to our previous study, 
Majorana fermion dark matter interacting with standard model fermions via a scalar mediator is the 
most favoured explanation of the galactic center excess when characterised by Bayesian evidence. 
We locate the most plausible parameter regions of this theoretical hypothesis using experimental 
data on the dark matter abundance and direct detection interactions, the gamma ray flux from the 
Galactic center, near Earth positron and anti-proton fluxes, the Cosmic Microwave Background, 
and galactic radio emission. 

We find that the radio data excludes the model if we include synchrotron radiation as the only 
energy loss channel. Since it was shown that inclusion of other types of energy losses lifts this 
exclusion we discard the single radio data point from our combined likelihood [34]. The rest of 
the data prefers a dark matter (mediator) mass in the 10-100 (3-1000) GeV region and weakly 
correlated couplings to bottom quarks and tau leptons with values of 10“^-1 at the 68% credibility 
level. 
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Appendix A: Bayesian Inference 


In this section we summarize the statistical background of our analysis. Let P{A\I) and P{B\I) 
denote the plausibility of two non-exclusive propositions, A and B, in light of some prior infor¬ 
mation, I. The probability that both A and B are correct is given by the conditional expression 


P(ABjl) = P(AjBI)P(BjI). 


(Al) 


Bayes theorem follows from the symmetry of the conditional probability under the exchange of A 
and B: 


P(AjBI) = 


P(BjAI)P(AjI) 


(A2) 


P(BII) ■ 

In this context P{A\I) is typically called the prior probability and represents the plausibility of our 
hypothesis given information prior the observation B. The likelihood function P{B\AI) indicates 
how accurately the hypothesis can replicate the data. The posterior probability P{A\BI ) quantifies 
the plausibility of the hypothesis A given the data B. The evidence P(Bjl) serves to normalize 
the posterior. 

For theoretical models with a continuous parameter 0 Bayes’ theorem can be recast in the form 

£(Bl0,I)7r(0,I) 


P(0jB,I) = 


(A3) 


e(B,I) ■ 

The posterior distribution can be used to estimate the most likely region of 0. The evidence is 
calculated via an integral over the full parameter space 


e(B,I) = f/:(Bj0,I)7r(0,I)d0. (A4) 

Je 

For more than one continuous parameters, 9i, marginalization is performed by integrating the 
posterior over various parameters in the higher dimensional parameter space 


m) 



(A5) 
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